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1 Total force from the electromagnetic field of the source 
1.1 Electric force 

The electric force on a dipolar particle of polarizability a e is is [?] 



F t e (r) = ^n{(a e E*(r)d i E j (r))}. 
The electric field from the source is [2] 

Er (r) = MoA^ 2 J G^ P (r, r', «) P 3 (r', w ) d 3 r 



(1) 



(2) 



where Gf,' P (r, r', w) is the electric Green function from the polarization currents and contains the trans- 



mission Fresnel coefficients {t s , p )- Gfj' P (r, r',w) is expressed on using the angular spectrum of plane waves 



i7i2-i7 2 z 



(3) 



where G y ' (K) = i (s^|i% s = Kxz,p i = - =p ^zj /(n»fco) , 7» = v^M^o - 

and rij = y/£ifM ■ ri\ (n-i) is the refractive index of the medium placed at z < (z > 0). 

The self-correlation function (E* (r) £Jj (r)) at the position of the particle, i.e., at r = ro will be 



Vi,V 2 



(E*(r )E j (r ))=M^ I G% p * (r,r',a,) G»f (r,r', u) wjp (t> v t> 2 ,u>) dV^r 



(4) 



W^i (r^r^w) — (P£ (r[) Pi (r 2 )) being the cross-spectral density tensor of the source polarization. Sub- 
stituting W^i (r'^r^u;) into Eq. (4) we have (E* (r) Ej (r)) explicitly expressed in terms of the angular 
spectrum 



(£*(ro)i^(ro)) 



9 2 4 



Jv 1 ,v 2 J Ki ,k 2 (27r) (2-7r) 



x e 



(5) 



We shall assume the correlation with a Gaussian profile [I] 

W<f> (riy 3>W ) = 5^ Hexp(-(|ri - r 2 |) 2 /2a 2 ) %/ (2^) 3 / 2 a 3 , 



(6) 



where 5^ (u>) is the normalized spectrum of the source and a represents its correlation length of the source. 
With a suitable change of variables, the integration over the lateral space coordinates gives a two-dimensional 
delta function S 1 - 2 ^ [Ki — Ka] which will simplify one of the integrations over K. The rest of integrals lead to 
a more simplified equation: 



(^M^(ro)) = \ /6 S< p n<<>) 



2, .2, ,4 
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e -2, oS7l ^ e -I ff 2 5R 72 2 d 2 K 
2^72 



l|2 n 2 fc 2 



7^2 > 1 



71 1 +K- 



t 2 a, 



72| 2 + X 2 



(7) 



Next, we will divide the force into conservative and non-conservative components. 
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1.1.1 Conservative force. 

The conservative electric force p ( r> cons — Keadi (E* (r) Ej (r)) /4 due to the electric field is given by 

T^e—cons £o£l < 



K = + x 



K=k 



»W T ,Ji/« * (P) H 

t s 21 r + \t 2 



4 ' — ' 4 (27r )V8 



x e 



-2z 97l . 



2^72 



| ,2' 
|72| 



(S) 



where only the third Cartesian component and the contribution is solely due to the evanescent modes. This 
integration is numerically resolved . 



1.1.2 Non conservative force. 

The non-conservative force F?' nc — Imadm(£<* (r) diEj (r)) /2 due to the electric field , is determined in a 
similar way, now calculating (E* (r) diEj (r)) instead of (r) diEj (r)). Hence, 



\l/5 v 



(27T) 



AT = fen 



K=0 
1 

23^ 



l"2 



|2 , up ,2 1 f,_ ,2 , T ^2 



|iiir + i4i 



(, 2 n,Q 



(|72 



(9) 



where only the homogeneous waves give a non-zero value which is constant for any value of r. This integration 
is resolved numerically. 



1.2 Magnetic force 

The magnetic force for a magnetodielectric particle is [lj 

Fr M = ^^{{a m H* (r) d^, (r)> } . (10) 
The magnetic field emitted by the source is Maxwell's equations, 

fff c (r) = -^/ Gf/(r,r', W )P fc (r',u;) < iV, (11) 

where the magnetic Green's function is 

Gf/(r,r» = ^/^(K)e*M')e^' ', (12) 

and (K) = — (^21^ — Sfc^ii^) - res ^ °^ ^he calculation is similar to the one described in Section 
1.1. 72 



1.3 Interaction force 

The force from the interference between the electric and magnetic dipoles is [1J 

F,r m (r) = - £o ei-^{(a:a m ) (E* x H) .} , (13) 

iZTT 

where Z = \J noH\/{eo£\)- Once we have defined the electric and magnetic fields emitted by the source [cf. 
Eqs. (2-3) and (11-12)], one can calculate the cross product (E* x H), 
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2 Total force from the electromagnetic field emitted by the particle 
induced electric and magnetic dipoles 

2.1 Electric force 

The field emitted by the electric (magnetic) dipole p (m) , after reflections on the source surface z = is 

Ef(r) = KiMxJ 2 f Gf j *(T,r',u)p j (i>,u)8(T'-T )<£ s T>, 

= fi fi 2 u 2 Gff p (r,r , u))pj (r ,w) , (14) 

= Zoiu> G H m ^ T , u) (w) ^ (15) 

c 

where the electric Green's function contains the properties of the source through the reflection Fresnel coef- 
ficients (r SjP ) 

G^(r,r» = ;|0Gt ? ( K )^" R ^ ni(z+Z,) . 

"'""♦(r.r'.w) - ^ /'^ G ^«- (K)el K(R-R') e ,7 1 (^'). (17) 
2 j (27r) 



G 



G^' P (K) = ^ (s^f 2 % +Pu r i2Piji) andGy ,mo (K) = i (pjrf 2 s,- - Sirf 2 p y ,). The superscript ° denotes 
that the electric field generated by the magnetic dipole has the same Green's function as the magnetic field 
radiated by the magnetic dipole with the interchange r s -H- r p . Notice also that in the Green function 
described above there is no free-space term. This is due to the multiple scattering of the dipole field with 
the source surface, (see for example [HI E| ) - 

The correlation function (E* (r) Ei (r)) at the position of the particle is obtained on considering that the 
cross-correlation between the electric and magnetic dipoles is zero, i.e., (p*mj) — 0, therefore 

(E* (r„) Ei (r )) = < (Ef* (r ) + E™* (r )) (Ef (r ) + E? (r„))) = (Ef (r ) Ef (r )) + (E™* (r ) E? (r )) , 

(18) 

each ensemble average being 



{Ef (ro)^(r )) = ^Mi^T (r 0j r 0|W ) G^ (r.r'.w) (p* k :(r ,cj) Pl (r ,cj)) , (19) 
(Ef* (r ) (r )) = Gf™(ro,ro^)G^(r,r', W )(™nro, W )m ; (r ^)), (20) 



and 



Pk(r,u>) = e Q eia e {uj)E l k nc (r ,n,w) , (21) 
p*(r,w) = a m (w)fff c (r ,r 1 ,w). (22) 

The correlation tensor (E™ c * (ro, i"i, ui) E\ nc (ro, ri, w)) has been calculated in the Section 1.1, where ri is a 
point of the half-space ocupied by source. After a somewhat protracted calculation one obtains that at the 
particle position one has that (Ef* (r ) Ef (r )) = 0, thus the only term remaining different from zero is 
(Ef (r )^(r )): 



(Ef (r )£f (r )) = /ig/*?w 4 |e | J | £l | 2 \a e \ 2 Gff (r , r , w) G^ p (r , r , w ) <£f c * (r , ri, w) £f c (r , r' lf w )> 

1 

(2^ 



2 2 4 I 1 2 I 1 2 I i2 Mo^ 1 c(-P) / \ 

A*oMi w M M |a e | 4 , ,i/ 5 S M 
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(23) 

After completing the whole procedure, we derive the conservative plus the non-conservative parts of the total 
electric force. 

2.2 Magnetic force 

The magnetic field emerged from either the electric and the magnetic dipoles is calculated by using the 
following substitutions [7]: 

E?{t) = -^tff(r), (24) 

= T^Ef(r), (25) 

p -> m, (26) 
r s O r p < (27) 

The procedure to obtain the force follows the same steps as described in Section 2.1 

2.3 Interaction force 

Once whe have characterized the electric and magnetic fields emerging from theelectric and magnetic dipoles 
[cf. Eqs. (14-17)], one calculates the cross product (E* x H) in a way similarto that of see section 2.1. 
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